Glass phase of two-dimensional triangular elastic lattices with disorder 
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We study two dimensional triangular elastic lattices in a background of point disorder, excluding 
dislocations (tethered network). Using both (replica symmetric) static and (equilibrium) dynamic 
renormalization group for the corresponding N — 2 component model, we find a transition to a glass 
phase for T < T g , described by a plane of perturbative fixed points. The growth of displacements 
is found to be asymptotically isotropic with u% ~ u\ ~ A\ In 2 r, with universal subdominant 
anisotropy u% — u\ ~ Ai In r. where Ai and A2 depend continuously on temperature and the 
Poisson ratio a. We also obtain the continuously varying dynamical exponent z. For the Cardy- 
Ostlund N = 1 model, a particular case of the above model, we point out a discrepancy in the value 
of A\ with other published results in the litterature. We find that our result reconciles the order of 
magnitude of the RG predictions with the most recent numerical simulations. 
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I. INTRODUCTION 



There is a lot of current interest in the problem of elas- 
tic lattices in the presence of substrate disorder. This is 
important in various physical systems such,as flux lattices 
in superconductopp, charge density _wavesa, colloids and 
magnetic bubblesQ, Wigner crystals!!, interface-.roughen- 
ing with disorder and solid friction of surfacesQ. An im- 
portant distinction is whether topological defects such as 
dislocations are present in the system or not. These de- 
fects may appear spontaneously because of disorder or 
not, or may be explictly excluded from the model (e.g 
flux lines in d=l+l geometryfftl, tethered networks with 
permanent bonds). Systems which are both elastic and 
periodic and do not contain topological defects are be- 
lieved to form glass states with many tnetastable states, 
diverging barriers and slow dynamicsBil3. There are very 
few analytical methods to study these states and their 
physics is far from being completely elucidated, despite 
help from extensive numerical simulations. When topo- 
logical defects are present the problem is even more dim- 
cult and very little is known. A crucial question of course 
is whether these topological defects, when allowed in the 
model, will appear spontaneously at low temperature be- 
cause of disorder. 

The general problem of a lattice on a . diso rdered sub- 
strate was addressed in a previous work£illi3 and it was 
predicted that in d = 3 and for weak pointlike (uncorre- 
cted) disorder a thermodynamic phase free of unbound 
dislocations should exist and be stable. In this phase, 
called the Bragg glass, translational order decays very 
slowly beyond a translational quastarder length R a (at 
most algebraically). These resultaL3lL3 being obtained 
by mapping the lattice problem onto a multicomponent 



version of a random field XY model the very same predic- 
tions (i.e absence of vortices and quasi long range order) 
obviously hold for the usual random field XY model as 
well. These results are at variance from .previous stud- 
ies which either argued forE-3 or assumeduJ the sponta- 
neoua|.genepation of dislocations, or did not address the 
issueEMO'El. Predicted consequences^ for the phase di- 
agram of type II superconductors, i.e existence of a low 
temperature low field Bragg glass phase of vortices under- 
going a transition into an amorphous state upon raising 
the field seem to be in agreementlia with the most recent 
experiments. The existence of a d — 3 defect free, Bragg 
glass type phase was confirmed in several numerical sim- 
ulations either in the contexi_pf superconductors^ or for 
the random field XY modelEj. Though a rigorous proof 
is still lacking there is also at present further theoretical 
support fojtJ.he stability of the Bragg glass: analytical 
variationaLrtJ and RG studiesu in a layered geometry and, 
very recently, a proof using improved scaling and energy 
arguments^. 

The Bragg glass is thus an example of an elastic glass 
phase with internal periodicity. In the physics of this type 
of glasses two dimensions play a particular role. Is-ad-i 
dition, further analytical methods become availableE3~E-3 
in d = 2. When dislocations are excluded by hand, it is 
believed that the lower critical dimension of these elastic 
glasses with internal periodicity is di c — 2, with a glass 
phase existing for T < T g in d — 2. For the N = 1 compo- 
nent model (random field XY model) thisas the result of 
the Cardy Ostlvmd (CO) RG calculation^. By analogy 
it can be arguedtj that the same holds in the case of the 
triangular lattice, but it has not yet been directly verified. 
A triangular lattice requires a fully coupled N = 2 com- 
ponent model. The CO N = 1 calculation could in prin- 
ciple describe a N = 2 decoupled square lattice, except 
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that such a lattice has usually a more complicated elas- 
ticity tensor which again couples the components. When 
topological defects are allowed, it was shown inE3 that for 
the N = 1 component model they are perturbatively rel- 
evant near the glass transition T g . As argued inE3 T g for 
the triangular lattice is well above the KTNHY melting 
temperature T m and dislocations should then be relevant 
near T g for the N — 2 triangular lattice as well. At low 
temperature however, much less is known about the im- 
portance of dislocations. The common beliefEl, which is 
by no means rigorously established, is that if dislocations 
are allowed, no glass phase will exist isurf = 2. This is 
also hinted at by the results obtained int3 on the simpler 
case of the random phase shift model, relevant tojdescribe 
a lattice with internal (i.e structural) disordered (a . sub -, 
set of the present problem). In these early studiesc3'E£l 
the high temperature phase with unbound dislocations 
was found to be rentrant at low temperatures suggesting 
the importance of topological defects— at low tempera- 
ture. However, it was pointed out int3, from a careful 
study of the CO RG flow, that at low temperature the 
scale at which the lattice is effectively dislocation-free 
(i.e the distance between unpaired dislocations) can be 
much larger than the translational length R a . Thus even 
in d = 2 the Bragg glass fixed point may be useful to 
describe the physics, as a very long crossover or maybe 
directly at T = 0. Furthermore, it was also pointed out 
ir£3 that the conventional CO RG will not be adequate at 
low temperature since it assumes a thermalized descrip- 
tion of the vortices and neglects important effect such 
as the pinning of dislocations by disorder (i.e that the 
fugacity depends on the position). A .similar idea was re- 
cently proposed by Nattermann et al.ES who reconsidered 
the simpler jeaadom phase shift model. They explicitly 
showed thatcJO was incorrect at low temperature and 
proposed a modified coulomb gas RG. It is still an open 
question how these new set of ideas and techniques apply 
to the present problem. 

In the present paper we study the statics and the dy- 
namics of two dimensional isotropic lattices interacting 
with point like disorder excluding dislocations. We use 
a (replica symmetric) Renormalisation Group (RG) ap- 
proach as well as equilibrium dynamics RG. Though we 
will not attempt to include dislocations, it is still interest- 
ing for several reasons. First, it is the natural continua- 
tion to d = 2 of the non trivial fixed point which describes 
the Bragg glass phase for d = 3. The only other analyti- 
cal approaches available to describe th^JSpagg glass phase 
are the replica-,variational methodOESEJ pid-. Fisher's 
functional RGa in a d = 4 - e expansiorJiilM Thus, 
the study of the present fixed point in d = 2 is also use- 
ful by comparison. Second, it is of general interest for 
systematic study of the d = 2 glass phases and allows 
to show that the lower critical dimension is d = 2. It is 
also a first step towards the study of the more difficult 
problem of lattices with disorder in d = 2 with disloca- 
tions allowed. Also, since it has more parameters that 
the N = 1 component model and can be similarly stud- 



ied numerically, it may lead to further useful numerical 
checks of the replica symmetric RG in this problem as 
well as the dynamics (see discussion below). Finally, it 
may be possible to realize such tethered networks exper- 
imentally. It was argued recentlyEl that internal disorder 
which breaks the internal periodicity of the network (and 
may occur in e.g polymerization) drive the system to an- 
other universality class. Even if this claim, which we 
will not investigate here, is correct, crossover lengths will 
probably be large for weak internal disorder. 

The N = 1 Cardy-Ostlund model, which is a subcase of 
our present calculation, has also generated a lot of atten- 
tion recently, for different reasons. The gaussian replica 
variational method (GVM) applied to this model, leads 
to a one step replica symmetry breaking (RSB) solution 
below T g and to mean squared displacements which grow 
as u 2 r~j Tglnr, a different result from the replica sym- 
metric (RS) RG prediction u 2 ~ A\ In 2 r. Note that the 
GVM, being by construction an approximation which ne- 
glect some non linearities has no reason to yield the Q*&pt 
result (see howeverEil). However it was also showr£3'L2l 
that the Cardy Ostlund RS-RG flow is unstable to an 
infinitesimal RSB perturbation at and below T g . The is- 
sue was thus raiseda of whether the RS-RG may miss 
some of the physics related to RSB. It is thus important 
to perform numerical simulations and carefully check the 
predictions of the RS-RG, as well as carefully define the 
observables to be computed. The numerical studies on 
the random phase sine-gordon model presently available 
show discrepancies, and their analysis is not satisfactory. 
IrO no effect of disorder was found, while inLJ the re- 
sults seem close to the predictions of the GVM. In the 
more recentEjO the results seem to agree qualitatively 
with the RS-RG but with an amplitude /rmich smaller 
than the the prediction for Ai quoted inE3E3 to which 
it was compared . Let us also point out the two very 
recent studieseElo of the random substrate SOS model 
(believed to be in the same universality class) directly at 
zero temperature which also yield u 2 ~ A\ In r. 

In view of the above interest, it is thus important to 
compute accurately the amplitude A\ predicted by the 
RS-RG. This is one of the result of the present paper. 
Our result, smaller by a factor of 4 from the previous 
resultE9'E3 isjdjffexent from any of the previously pub- 
lished resultsHEjo and, we argue, is correct. This seems 
to reconcile, at least in order of magnitude, the result of 
tlio-RS-RG with the result of the numerical simulations 
o£jO. A fiuejj-and more precise analysis than the one 
performed inr^r 7 ! is however needed, including predicted 
finite size scaling corrections, since a more careful treat- 
ment of the effects of RSB may— reveal that deviations 
from the RS-RG pesult are smallc 2 ] e.g only in the ampli- 
tude of the log 2 rc3. Another consequence of our results 
is that the distribution of rescaled displacements w/lnr 
should be gaussian at large scale, which we propose as a 
useful numerical check of the RS-RG method. 

Along these lines, one notes that the equilibrium dy- 
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namical RG studied here, because it obeys by construc- 
tion the Fluctuation Dissipation Theorem (FDT), gives, 
as we illustrate explicitly, identical answers for static 
quantities than the replica symmetric static RG. (it thus 
provides for us a useful check on our statics calculations). 
As discussed above, it is possible that the static RG does 
need a RSB structure (either strong, i.e in the ground 
state as in mean field - or weak, i.e only in the low ly- 
ing excitations) and a proposal for its construction was 
made inE3. Also, we interpret the recent resultsEj on the 
simpler random phase shift model, i.e its mapping on 
some version of Dcrrida's random energy model, as an- 
other hint that RSB may be important in these models. 
If this is the case, it is then obvious that an equivalent 
statement should exist in the dynamics, i.e that one may 
need to construct a RG which violates FDT and replaces 
it with a generalized FDT structure, a progra m .. W |hich 
was successfully carried in mean field dynamicsoc3 We 
will not address this question further here but will lucidly 
concentrate on the predictions of the replica symmetric 
RG which is certainly a good starting point in this prob- 
lem. 

The paper is organized as follows. In Section (|J) 
we define the model for a triangular lattice on a disor- 
dered substrate. In Section (III) we study the statics 
of the problem. First we use a mapping onto the repli- 
cated Coulomb gas and derive the RG equations (Section 
( III A ). Then we anal yze th e RG flow and compute the 
correlation functions (IIIB). The res ults fo r the tr ian- 
gular lattice can be found in ( III B 1 ). In ( III B 2| ) we 
give the results for the N = 1 Cardy Ostlund model 
and review the discrepancy with other publish ed res ults. 
In Section (IV) we study the dynamics. In ( IV A ) we 
carry first order perturbation theory and in (IVB) we 
obtain the results for the dynamical RG. We then com- 
pute the dynamical exponent, and obtain the scaling 
behaviour of transport quantities at the transition. In 
particular we identify an "effective critical force" at the 
transition. Section (|y|) gives the conclusion. Technical- 
ities are releguated in the appendices. In Appendix (^) 
we discuss in detail the regularization, obtain the RG 
equations and show explicitely the cutoff independence 
of some universal ratios. In Appendix ([b]) we perform 
the RG directly on the static replica effective action. In 
Appendix (III) we perform the dynamical RG to second 
order. In Appendix (IV) we study some consequences of 
the statistical tilt symmetry. 



positions i?° to new positions Ri = P® + u(R®). Further- 
more the connectivity of the lattice is fixed: each site 
will have exactly six neighbors (see fig 1). This can be 
realized in principle by considering a network of identi- 
cal and permanent nearest neighbor bonds (tethers) with 
the topology of a perfect triangular lattice. Dislocations 
are thus excluded in this model by construction. As dis- 
cussed in the Introduction, the model can also be rele- 
vant at length scales and in regions of the phase diagram 
where dislocations can be neglected. 




FIG. 1. Representation of the elastic lattice and its recip- 
rocal lattice vectors of minimal modulus. 

The interaction energy can be described in terms of 
the local displacement field u(f) by the harmonic hamil- 
tonian 



H 



d 2 q 



(1) 



where §ij(q) is the elastic matrix of the 2D lattice. 
In this paper we will consider isotropic (i.e triangular) 
lattices which can be described by only two independant 
clastic coefficients : 



$iM = c 11 q 2 P l L i 



(2) 



where P?j = q t qj = qiq d jq 2 and Pj- = S, tj - q^ = 
(q±.)i{q±)j are respectively the longitudinal and trans- 
verse projectors. 

The interaction of the lattice with the impurity disor- 
der of the substrate is modelled by adding to the hamil- 
tonian: 



H v = d 2 x p{x)V{x) 



(3) 



II. THE MODEL 



We consider a set of identical "atoms". The interac- 
tions between them is such that the ground state config- 
uration in the absence of disorder is a perfect triangular 
lattice (see fig. 1) of equilibrium positions and lat- 
tice spacing ao- Both thermal fluctuations and disorder 
will result in displacements of the atoms from their ideal 



where p(x) is the local density of atoms, p(x) — 
^ 2 \x— Ri). The random potential is gaussian with 

short range correlator V(x)V(x') = h(x — x'). The sym- 
bol — denotes average over disorder, while (..) denotes 
thermal averages. 

-This model can be transformed, as described in details 
ir£3 (see Section II-B), into the following random phase 
model: 
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z = 



J d[u]exp - i (H + H dls ) 



T 



Hdis 



*2 



aij(x)(diUj + djUi) 



(4) 
(5) 



Coulomb gas formulation. In the second part we compute 
the static displacement correlation functions. 



A. derivation of the RG equations from the 
Coulomb gas 



The first terms corresponds to random localpStresses 
and comes from the q ~ part of the xiteordeniH. The 
effect of only this term was studied ino (in the pres- 
ence of dislocations), and would also arise in a problem 
of structural disorder (with no substrate). The second 
term originates from the Fourier components of the sub- 
strate potential Vk close to one of the reciprocal lattice 
vectorsE3. The random field is a phase distributed uni- 
formly over [0, 2tt] and satisfies: 



( e *W*(*)-*„'(*')) = S v y5W(x - x') 



(6) 



Note that we are keeping only the terms which are rele- 
vant in the RG sense near T g . There are additional higher 
harmonics terms which are irrelevant in d = 2 near T g . _ 
There are also higher non linear terms which are small in 
the elastic limit \u(Ri+{) — u(Ri)\ <C a where the deriva- 
tion of (|]J) is valid and also irrelevant by power counting. 
These term will cotrect the bare values of the relevant 
coupling constantsEa 

The random stress field has the general correlator: 



o-ij(x)a k l(x') = — [(An 



(7) 



+A 66 (d ik dji + SiiS jk )] 

The bare value is A 6 6 = and An - P 2 h k=0 /T (fromB 
) but they will flow under renormalization. 

We will study both the statics and the dynamics of 
this model. The statics of this model can be studied 
ing replicas to average over the disorder potential V(x) 
The replicated model reads: 



Z n = 



1 



d[u a )exp --[m 



(n) 



H, 



H, 



(n) 



-(AuP, 



■ A^)] q 2 u1{q)u){- q ) 



(8) 
(9) 



El 

T 



-5 E E cos ^ 



T{x) -u\x)} (10) 



g is related to the fourier coefficient A_r- and the con- 
stant density p : g = plh Ko /T 2 . We will also use the 
dimensionless coupling constant g — ga 2 . 



III. STUDY OF THE STATICS USING A 
COULOMB GAS FORMULATION 

In this Section we study the statics of the model. We 
first derive the renormalization group equations from a 



The model (||) is the analogue for a two-component 
field of the two dimensional random field XY model 
whose RG equations have been derived using the 
Coulomb Gas approach by Cardy and OstlundEIL Jiut 
contrarily to them, our Coulomb gas formulatiosc3E£l is 
obtained by first using a Villain approximation!!!! before 
averaging over the disorder. This results in a natural 
coupling between the continuous replicated displacement 
field u a (x) and vector charges n®(x) = n a (x, v).K v where 
K v is one of the three (f = 1, 2, 3) reciprocal lattice vec- 
tor of minimal modulus (see fig 1) and n a (x,v) 6 Z™ 



Z» = J d[u a ] ex P \-f H o-i / 

[»?(*)] v ' ' 



nl.u a (r) (11) 



with the condition \/{x,v) '^ a n a (x,v) = 0. The RG 
equations will be derived to lowest order in the charge fu- 
gacity y/g, since higher order operators are irrelevant near 
the critical pointH : thus we can only consider charges 
of the form n a (x, v) = <5° — 8% = 8® * where a and j3 
are replica integers between and n. The effective fu- 
gacity of these minimal charges is g. Integrating over the 
smooth field u(x) one recover a 2D vector Coulomb Gas 
with charges having both a spatial and a replica struc- 
ture, and whose action is (omitting the fugacity term): 

z v,v> 

T 



E / . VAWtfi?- f)ri b vd (f) (13) 



where the interaction V(r) is obtained by a Fourier 
transform. We have incorporated the replica off diagonal 
terms in . 

In this section, as is usualy done, we will-useE3, instead 
of the full interaction, its asymptotic formula (r ^> a) (see 
Appendix A) : 



a 



1 

2"* 



^ (14) 



, ab 



ab 



and ref = (egg 1 + c 11 i )"" /4tt, n, 2 
We will use the notation nf 3 = K\8 a b 



, ab 



■66 



/An. 



Ki and (cn) a 
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We won't deal with possible appearence of dislocations 
in this model. This will lead-to an electromagnetic vector 
CG which neads more workE.il. This simplification allows 
us to derive RG eouations by working on the 2-point 
correlation function^ instead of the replicated partition 
functionE2rc3. These two schemes are equivalentE3 to or- 
der (at least) g 2 . The renormalised elastic matrix can be 
defined as 



, a 6 



.2 



-1 
66, R 



\ab pi 
) r i'i 



lim — 

9^0 T 



{u1{q)u b A-q)) 



(15) 
(16) 



To be consistent this definition must be independant of 
the direction q : this is simply a direct consequence of 
the isotropic nature of the 2nd order invariant tensors 
on the triangular lattice. The explicit expression for the 
renormalised elastic coefficients naturally follows : 



(c- 1 -c~ l \ 

l L ll,R c ll J 



-Tlim q 2 {ni :l {q)ni, (-q)) s[n] P, 



(17) 

bd 



where I, m are cartesian coordinates, and we used the 
(j^uj . The same holds for Cgg 1 by 



notation n c v l = n c 



simply replacing Ph by P/, 



(P.a 



(a,P) 



(a) 



(a,P)" * 
. R 

(b) 



(y.a) 



(a,p 
t5 



(P.Y) 

i 

(c) 



FIG. 2. Representation of the three different configu- 
rations that contribute to the correlation function of the 
Coulomb Gas charge n. Configuration (a) is of order g 2 
whereas (b) and (c) are of order g . 

Replacing <I> _1 by its expression (j|) and using the 
Fourier transform of the charge correlation, we obtain 



(18) 



This equation is the starting point of our RG study. 
Note that itis a definition of the renormalised elastic co- 
efficient exact to all order in g. We now proceed as follow 
: first we expand perturbatively in g the right-hand side 
of this equation. One can then apply a coarse-graining 
procedure which leaves the form of this perturbative ex- 
pansion unchanged, thus defining coarse-grained elastic 
coefficients. To order g 3 only three terms appear (see 
fig. 2) : (a) of order g 2 and (b) and (c). Corrections 
coming from the coarse-graining of configuration (a) de- 
fines the coarse-grained elastic constant cT/j 1 (Z) to order 
g 2 and g(l) to order g, where I is the scaling 



came- 

tcr. In the spirit of Coulomb Gas renormalisatior£3, this 
corresponds to the annihilation of charges (screening of 
the coulomb interaction), whereas the fusion of charges 
is given by the contributions of configurations (b) and 
(c) when \R\ <C \r\ or when \R — r \ <C \ f\. This gives the 
correction of g(l) to order g 2 . 

First, let us develop the two point correlation function 
(nl i (0)n^ / j (r)) to order g 3 : it is given by the three 
contributions (a), (b), and (c) (see fig 2) : 

«z(0)r?^ m (r*)) sw = g 2 A 1 (6, r) + g 3 A 2 (6, r) (19) 

where the first term corresponds to the only neutral 2- 
charges configuration : 



-S a 



(20) 



and configurations (b) and (c) contribute to the factor 
A 2 : 



« n 

A 2 (0,?) =2/ EE 6 c a , f3 6 d ajf3 (K„) l (K„>) m e- S <> 

• )R v=£v' 
« n 



Sc 



In these expressions, S a ,Sb,S c are the respec- 
tive actions of the three configurations : S a — 

2T\K\ 2 (k,! \n{r/a) - n 2 {K.f) 2 + k 2 /2) . The two others 

can be easily expressed in terms of S a . We can put this 
back in the expression (|l8|) of c^ R . The angular inte- 
gral of the g 2 term can be performed, using the parameter 
a = n 2 \K\ 2 T : 

, Sir 



K 2 —(±h(a)+2I (a)) 



(22) 



where Iq and I\ are the modified Bessel functions. The 
function B\\ appears in the expression of c^ R while Cgg 1 R 
involves Bqq . Using the following results on the replicated 
charges 
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^5X = -2(n^-l) 

8 a a An = ( 2 - n ) ( n6ab ~ !) 



where the fonction B g (a) is given by 



B 9 (a) = / di?e- Q / 2 [(2-n)e 



(if, A) 



(27) 



-2a cos(K,_R) cos(-RT',il)l 



li 



(23) 

3-2|7?| 2 KiT 



and the definition g — ga 2 , we can express c lx fl as 

Ki 1 ) (0 

l)Bn(«) f-(-) 

c ll «'|» s l>o a 

The integrand in the last term is given by 

d 2 R3 /r\2 



= 27r[(2-n)/ (a/2)-/ (a)] 
Integration over r leads to 



-K 2 SA(r) = 2B g (a)B 11 (a) 



a Va/ 



This is exactly the expected form, which can be written 



as a correction to g 



where 



A(r) 



\R\>a 



a 2 8 



(-) (cos2(^,f) + 2) 



d(g 2 ) = -dl g 6 2B g {a 



B g (a) = -2tt (I (a) + (n - 2)/ (a/2)) 



(28) 



x ((2 - n)e" 



We can then apply the coarse-graining procedure to 
this self-consistent equation. We thus rescale the hard- 



The final equations, obtained after the n — > limit, are 
d en d c 66 



core cut-off from b to a 



adl. The first 



integral split into two terms (r) = \K\ 2 niT) 
dr /r\3-2 t; 



d I dl ° 
= g 2 TB n (a) 



d An(Z) _ _ 2 



d / 
dA 66 (Z) 



f — (L\ 

J a a ^ 



a Va 



dr / r\ 
i a Va/ 
dr /r\ 3 - 2? 7 
a Va/ 



dZ + e'( 4 - 2 ") 



3— 277 



d / 



<7 TB 66 (a) 



^| = (2-^ 2 / tl T)s-ff 2 B 9 (a) 



(29a) 
(29b) 
(29c) 
(29d) 



Hence the first terms contribute to c n (I) while the 
second renormalises g(l) : 



d [cn ] (0 = ~g (nS ab - l)TBu{a)c^ 2 dl (24) 



d 5 2 (0 = (4 - 2\K\ 2 KiT)g 2 dl 



(25) 



For the purpose of this paper, we are only interested 
in the contribution of the last term of eq. ( |23|) to dg 2 . It 
corresponds to configuration (b) and (c) of (fig. 2), when 

a < \R\ < a(l + dl) or a < \R- r\ < o(l + dl). 

The renormalisability of the model that one can di- 
rectly see in our procedure ensures that this correction 
can be written as a a contribution to g 2 (l) of order g 3 . 
First of all we check these renormalisation properties of 
the hamiltonian : 

lim e-^^) =e~ SJf) e -i<x(Ki.R)(K2-R)-% 

|i?|->a 



with g = ga 2 . From now on B g (a) denotes its value at 
n = 0. 

It is useful to compare these equations with the one 
for the n = 1 component model of Cardy and Ostlund. 
To obtain them we set en = egg = c (i.e a = 0), 
An = A66 = A and consider the two reciprocal lat- 
tice vectors (instead of three) of a square lattice in the 
interaction term with K 2 = 1 . This gives two decoupled 
n = 1 component model with RG equations obtained by 
Bn = Bqq — ► 7r and drop the term not proportional to 
n — 2 in the g 2 correction. 



dc 
~dl 
dA(l) 
dl 







d I ZTTC 

the transition temperature is thus = Anc. 



(30a) 
(30b) 
(30c) 



lim e - s c( F > R ) = e ~ Sa ^ e"*^^) 2 -! B. Analysis of the RG flow and static correlation 

>o functions 

The corresponding correction is TT . , n , ,. , 

We now study the KG now and compute the correla- 

3 /r-\ 2 / „ s tipHfc-iIn the case of the N = 1 l»£jdel this was first done 

SA ^ = l\a) ( cos2 (^) + 2 J e aB M) (26) irSil and later reconsidered i: " 



G 



1. triangular lattice 

The flow defined by the above RG equations (^) has 
similarities— with the one obtained for the random field 
XY modeM There is a transition at T g = 2/(K 2 n 1 ) = 
8irK~ 2 ciic 66 /(cu +cqq) and one has K 2 = l&ir 2 /ia 2 for 
the triangular lattice. 

In the high temperature phase T > T g the disor- 
der renormalizes to zero. At low temperature T < T g 
the coupling constant g(l) converges towards a pertur- 
bative fixed point g*. Introducing the reduced temper- 
ature r = (T g - T)/T g , one has: g* = 2T/B g (a) with 



B g (a) = 27r(2/o(a/2)-/ (a 
on the value of a, defined in 
ate at T g since we work near T g 



It depends continuously 
and which we can evalu- 
Using the above value for 
T g one finds that a(T g ) = a g = 2(cn — C66)/(cn +C6e) = 
2(1 + er)/(3 — a) where a is the Poisson ratio defined as 
usual as a = 1 — 2cee/cn = A/(A + 2/i). Thus we find 
that the obtained fixed point depends continuously on 
the Poisson ratio. Since the Poisson ratio is not renor- 
malized one has now a plane of fixed points, rather than 
a line in the case of the N — 1 Cardy Ostlund model, 
parametrized by the temperature r and the Poisson ra- 
tio a. 

One must check however that the perturbative fixed 
point does indeed exist, i.e that B g (a) > on the allowed 
domain of variations of a. This is indeed the case since 
one finds that B g (a) > as long as a < a* « 2.218. This 
condition is fullfilled since cqq > implies that a < 2. 

At this fixed point both An (7) and A 6 6(0 grow to non 
perturbative values. [-This is not a problem since due to 
statistical symmetryE3 they do not feedback in the RG 
equations. This produces however a change in the corre- 
lation functions (see e.g Appendix IV). 

We can integrate the above RG equations and the so- 
lution reads: 



For weak disorder g <C 2r/B g it corresponds to the 
length £ ~ R a ~ (2t/# ) 1/(2t) > a at which trans- 
lational order decays asymptotically (also equaLto the 
Larkin length R c for this lowest harmonic modelE, valid 
near T g ). Note that one can define another length scale 
£ c defined as the crossover length between a log(r) be- 
haviour of the mean squared displacements (see below) 
for r < £ c to a large r regime which corresponds to a 
log 2 (r) asymptotic correlation. This length diverges ex- 
ponentially as £ c ~ exp(^cst/r 2 ) at the transition r = 0. 

We now compute the correlation functionsEa. We will 
be interested in the following correlation functions: 



B«(r) = ((ui(r) - u*(0))K-(0 - ^(0))) (36) 
' hq (1 - e^iuMM-Q)) 

L i T3_(~\ r>T 



(2n) 2 

B L (r)P^ + B T (r)P i:j 



where the projectors are defined by = fjfj and 

Pij = (f s ±)i(f : ±)j- The longitudinal correlator is thus 
given by 



B L (r) = 2 



d 2 q 

WT 2 



(l-e^P^uMM-Q)) (37) 



The transverse correlation function follows by replacing 
P L by P T . 
We define 



P y (q,A 11 , 66 (0),g(0)) = (uM^(-q)) (38) 
Using usual dimensional scaling relations, we can write 
Tij(q, A Ul68 (0),ff(0)) = e 21 ^(e 1 ?, A Ui6a (0,ff(l)) 



9(1) = 



go e 



2tI 



1 + X(e 2rl - 1) 



(31) 



A n ,66(0 = A(0) + DnfiG I ^ ' 



2t 



ln(l+ X (e^-l)) (32) 



-x(i-x) 



2tZ 



1 + x(e 2Tl - 1) 



where we defined -Dii,66 = 4r 2 T 'i?n. 66(a) / B 2 (a) and 
X = g B g (a)/(2T). Note that X = So/0* for r > 0. Thus 
at large I one has in the glass phase (r > 0): 



We then choose the scaling parameter I such thato 
qe l = 1/a. The large r behaviour of the correlation func- 
tion thus corresponds to the limit e l — > oo. The RG flow 
approaches its fixed point in that limit : g — > g°° w g*. 
Assuming that r,y(l/a, A 1166 (i), g*) can be evaluated 
perturbatively in g* near T c We can then evaluate the 
correlation function : 



lW,Aii lM (0),fl(0)) = 

T / An(0 A660) 

„2 I ^ C H + „2 + l C 66 + 2 ^ij 

y \ c n L 66 



(39) 



,g(0-r + O(e- 2ri ) 



(33) 



(VP 



Aii, 66 (0 = Du,66 In ( y ) + A u ,66(0) + G(e- 2r ') (34) 



where we have introduced the length: 
'1 - x~ Inx' 



£ = a exp 



2r 



(35) 



The correlation function then takes the following form 



near T g : 



*o T 



(uMu^-q)) * --fln(«0 -^Ai + -2^^ (40) 

9 \ c ll c 66 / 

The angular integrals can be easily performed using 
the formula: 
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1 

2-k 



£ / dq(l e i nP°P§ = - (1 - J Q (qr) + E af} J 2 {qr)) 



(41) 



proportional to \n{r/a) (e.g a term like A(0) \n(r/a)). 
Therefore when comparing the RG predictions with e.g 
numerical results, one should be careful in identifying the 
different components of the spatial correlations. 



with a, (3 = L,T and ell = ett = 1 and elt = £tl 



-1. 

This gives: 

B L ,T{r) = ^{B lx +B m ) f ^ln(l/(«0)(l-Jb(«r)) (42) 
2?r J q 

±±(Bn-B u ) J ^ln(l/(qO)J 2 (qr) 

where -Bn,66 = ^g-Dii.ee/cn 66- The integrals can be 
evaluated as: 



600.0 



1 

2^ 



V^(Wo(*r)) 
o 9 



= -^ln 2 (I 



1 



• O ( In ( - 

- a 



2tt Jq q 



0(1) 



The difference between the expressions of Bl and £>t 
is the sign of the second order Bessel function ^(gr). 
Both these correlations have the same leading term in 
the limit of large r, given to the lowest order in r by: 



B L (r) ~ B T (r) 



(43) 



with a universal coefficient (the cut-off dependence 
drops out) 6(a) = (Bn + £?66)/(47r) Another univer- 
sal quantity (in the limit of large r) is the difference 
Bt{t) — Bl{t) since the leading integral, with all its cut- 
off dependent corrections in In r/a cancel exactly. One 
obtains: 



B T {r)-B L (r) 



b(a)_ r2]n/ r 



K 2 



(44) 



The coefficient of the log is a universal function of Cn 
and C66, 6(a) = (B 6& — Bh)/(Att). Specifically one finds 
(see fig. 3 for a plot of these functions) : 



6(a) = 6 
6(a) = 6 



2J (a)(l + a-)-aJ 1 (a) 
(2/ ( a /2)-/ (a)) 2 

Ji(a)(l + ^)-2af (a) 
(2/ (a/2)-/ (a)) 2 



(45) 

(46) 
(47) 



where a = a g = 2(c n - c 66 )/(cn + c 66 ). 
Here it is interesting to remark that in addition to the 
two terms (|43|,[44|), Bl^t^) contains a non universal term 



400.0 



300.0 



400.0 



- 300.0 




600.0 



100.0 - 



100.0 



FIG. 3. Plot of the numerical coefficients of the correla- 
tion functions Bt,l(t) and Br(r) — Bl(j) : 6(a) and b(a) 
as function of the elastic ratio a — 2(1 + c)/(3 — a). We 
have also plotted the coefficient of the dynamical exponent : 
0(a) = (z — 2)/t. 6 co , the corresponding value for the planar 
random field XY model (the Cardy-Ostlund model) is shown 
as a reference. See the text for the details. 



As we discuss in the Appendix IV the distribution of 
rescaled displacements u/lnr should be gaussian. Thus 
we have also obtained the decay of the structure factor 
at large r: 

exp(-i6(a)r 2 ln 2 (^)) (48) 



(J K t r \ — e iK .u(r) e -iK .u(0) 



the corresponding result for the N = 1 CO model 
structure factor is given in |l0| of the Appendix |TV| . 
Though the In 2 r term is isotropic there is subdominant 
anisotropy. For instance one has for the decay in different 
directions: 



C'K (r || K ) 



f 6(q)t 



(49) 



C Ko (r± K ) 

which is the analogous for d — 2 of Eq. (4.32) of RefM 
High temperature phase 

Finally, the high temperature phase is characterised by 
g(l — » oo) = and a non universal value of Aj(Z) — ► A?°. 
From (|3l|) we find 
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^11,66 - ^11,66(0) 



Di 



1,66 



2r 



(ln(l -x) + X) (50) 



which depends on the initial bare values of A but remains 
finite at the transition r — > CP. 

Since the renormalised theory is gaussian, the corre- 
lation function is straightforwardly given in the limit of 
large r by: 



((u t (r) - Ul (0))(u 3 (f) ~ Uj (0))) 

—Sij In (- 
7r \a 



(51) 



A°° A°° 
fe 1 + ^) + + 3*) 



li 



-66 



2. Cardy-Ostlund model 

In the case of the Cardy Ostlund JV = 1 problem, de- 
fined in its replicated version, as: 



H 



E 



(yp( cS ab - A a6 )Vw Q Vu 6 - gcos(K(u a - u b )) 



one finds, using (BO 



B(r) = ((«(r) - u(0)) 2 



(53) 



where, we recall r = (T g — T)/T g . As announced in the 
Introduction, this result is different from all the results 
previously published (to our knowledge). Since it is im- 
portant for comparison with several existing numerical 
simulations, we now discuss in more details this discrep- 
ancy. 

Our result is smaller by a factor of 2 from the original 
result (5.24) of Ref.cJ. It-is smaller by a factor of 4 from 
the result quoted in RefjEa, and still by a factor of 2 from 
its later corrected valueO (it4s also larger by a factor of 
4 from the result obtained inED). 

As discussed above the RG equations which read at T g , 
dA/dl — cig 2 , dg/dl = c^g 1 are not-universal but the fol- 
lowing amplitude ratio, defined ir£3 R = T g c\/ {cc2) 2 ^s 
universal. We find here that R = n. (as well as in Ref.Q). 
This is also the valued w e inf erred from the static and 
dynamic RG equations oSEJ and we thus agree with 
their RG equations. This value was R — 2n and thus 
incorrect inE3 but later corrected back to R = ttzB. 

The origin of thft-discrepancy between ( |53"| ) and the 
result (5.24) of Ref.Ej thus lies in the calculation of the 
correlation function. It can probably be traced to the 
algebraic mistake between equations ife-23) and (5.24) of 
Rcf.a. Finally the discrepancy withED lies first in their 
RG equations with N = 1 (we extracted their amplitude 
ratio as being R = it/ 2) and later in a factor of two 
in the calculation of the correlations (correcting for the 
additional misprint between formulae (13) and (14)). 

As discussed in the Introduction this improves the 
comparison betweenjiumerical simulations and the RS- 
RG. Indeed in Ref.E3 a result smaller by a factor of 



5 thi 
thanE 81 

result; 



:8 I was found, and we find here a result smaller 
y a factor of 4. Concerning the recent quoted 



of the numerical determination of the ground 
state, we find that it is about twice the naive continua- 
tion of our amplitude to T — 0, i.c setting (t = 1) in the 
above result (p3|). 



IV. DYNAMICS 



In this section we study the dynamics of the triangular 
lattice on a disordered substrate. We study directly the 
dynamics of the original random sine gordon model (|^). 
The method we use is to perform the renormalisation on 
the dynamical effective action associated to (||) and com- 
pute dynamical quantities. We find that this method 
is more convenient, and leads to more tractable calcu- 
lations, than the methodEEl originally used by Gold- 
schmidt and Schaub to study the dynamics of the Cardy- 
(5§^tlund model. In order to obtain the dynamical expo- 
nent z one needs, in addition to purely dynamical quanti- 
ties, some information on the statics. It would be poten- 
tially erroneous to attempt to extract from the previous 
Section the necessary information on the statics, since 
we used in the previous Section a different method with 
a different regularisation scheme. Instead, the consis- 
tent approach is to use the same method and the same 
regularisation scheme in the statics and the dynamics. 
We will thus carefully reobtain the results for the stat- 
ics within the same dynamical RG. This is possible, as 
we will demonstrate, because we are studying the equi- 
librium dynamics. Using the fluctuation dissipation the- 
orem (FDT), which then holds exactly by assumption, 
one can then reobtain the statics. We will also show 
(see Appendix [b|) that the same effective action method, 
but applied to the replica symmetric theory gives the 
same results for the static as the FDT equilibrium dy- 
namics. Also clearly demonstrate how the cutoff depen- 
dence comes in the RG equations and how at T = T g 
the physical amplitudes become universal, as detailed in 
Appendix 



A. perturbation theory on the dynamical action 

The dynamics of the model (|^) can be described by a 
Langevin type equation: 



d 



5H 



Su a (x,t) 



+ C(x,t) 



(54) 



with (C a {x, t)^(x', t')) = 2r]T5 a!3 d(x - x')5(t - t') be- 
ing the thermal noise and rj is the friction coefficient. The 
hamiltonian H is the sum Ho + Hdis defined in (||) . The 
equation of motion reads, specifically: 







V ^u a = c 66 W 2 u a + ( Cll - c 66 )d a dp U P (55) 
fi( x ) = -ZTy/g^K? sin (K v .u{x) + <f> v (xj) (56) 



T 



(57) 



A convenient method to study the dynamics is to use 
the de Dominicis-Janssen (or Martin Siggia Rose) gen- 
erating functional. Using the Ito prescription it can be 
readily averaged over disorder. The disordered averaged 
functional reads: 



Z[h, h]= J DuDue~ s ^ +hu+ihil 
S[u,u] = So[u,u] + S 2 [u,u] + Si n t[u,u] 



(58) 



S [U,U}= I iu^ qt ( V d t + Cll q 2 P^ + c 6(i q 2 P^)u^ t (59) 

Jq,t 



rjT / (»<%)(«%) 

Jr,t 



S 2 [u,u] = -- / (iu« t )(iut \ t >)q 2 



(60) 



Sint [u, u] 



x (A 66 P^(q) + A lx P^(q)) 
{iu a rt )(iu p rt ,)A a P{u rt -u rt ,) (61) 



where the correlator of the pinning force is 
fi(r,u rt )fi(r',u r ' t ') = A af3 (u rt - it r 't')- Z,From this 
functional one can obtain the disorder averaged corre- 



lation function C„ 



,a,/3 
rt,r't' 



u" t u^, t ,} and response function 



8{u? t )/5hP, t , which measures the response to 



R 

a perturbation applied at a previous time. They are ob- 



tained from the above functional as C, 

? "/3 _ /„,a-~l3 



and R"f r , t , — (u" t iu^., t ,) s respectively. Causality im- 
poses that Rrt^r't' = for t' > t and the Ito prescription 
imposes that R r t.r't — 0. All correlations iiii vanish. 
We will assume here time and space translational invari- 
ance and denote indifferently C r t, r 't' = Cr—r',t—v an d 
Rrt,r't' — R r -r' ,t-t' by the same symbol, as well as their 
Fourier transforms when no confusion is possible. Note 
that in this problem C rj _t = C r _t- 

In the absence of disorder the action is simply 
quadratic S — So and the response function is thus (for 
t > 0): 

Rft = P^(3)/ie- Cll9V 0(i) + Pl ( M^ c ^ 2 ^9{t) 

(62) 

we have introduced the mobility /i = l/i]. The corre- 
lation function is: 



c; 



a/3 



c\\q z c m q* 



-c 8e q fit 



(63) 



They satisfy the fluctuation dissipation theorem 
(FDT), i.e that: 



R af3 - 



(64) 



In the presence of disorder one will study perturbation 
theory expanding in the interaction term Si n t using the 
quadratic part So + S 2 as the bare action. The disorder 
has a quadratic part S 2 which is purely static and is im- 
material in the perturbation theory. Indeed the response 
function of Sq + S 2 is identical to the one of Sq and the 



correlation function is changed as: C^f 
with: 



cf t +a 



a 



a/3 _ 



statqj 



T^PZM+T^-PtM 



a/3 
statqj 



(65) 



which is purely static and does not appear in any di- 
agram of perturbation theory. Thus for any practical 
purpose one can consider that So is used as the bare ac- 
tion. 

We will perform the calculations on the general form 
for the correlator of the disorder : 



e iK.; 



(66) 



A' 



In the case of the triangular lattice for later conve- 
nience the sum is over the six reciprocal lattice vectors 
K = 1, 2, 3, —1, —2, —3. We will specialize at the end to 
the case of interest here, i.e a triangular lattice with the 
model introduced in Section II: 



A af} {u) = 2gT 2 



K"=l,2,3 



K a Kp cos(K.u) 



(67) 



where we denote by 1,2,3 the three reciprocal lattice vec- 
tors K w with v = 1, 2, 3 of (§). 

In fact, the condition A^ = AnKaKp is a potential 
condition which is preserved under RG and describes the 
dynamics of the model studied in the previous Section. 
With no loss of generality A(u) is an even function and 
thus Ak — A-k and Ak is real. The model of Section 
II is thus obtained for Ak — gT 2 , and more generally 
the relation with the statics studied in Appendix C is 
A K = 2T 2 g K . 

To establish the dynamical RG equations we will com- 
pute the effective action in perturbation of Si n t using 
So + S 2 as the bare theory. The calculation to second or- 
der is performed in the appendix and we will only quote 
the results. Here we compute the effective action r[u,{t] 
to lowest order in the interacting part Si n t'. 

T[u, u] — Sq[u, u] + S 2 [u, u] (68) 
+ (S in t[u + 5u,u + 5u])suM + °( s int) 

where the averages over Su, 8u are performed with the 
action 5*0 + 52. The calculation gives: 
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I>,u] =S Q + S 2 - f R^ rt ,(iu? t )(A a ^(u rt - u rt ,)) T[u,u} = S + S 2 + [ (tt%)E(«rt - u rV ,t - t')u% 

Jrtt' Jrtt' 

11 (iu^)(iu rt ,){A^(u rt -u rt ,)) (69) ~ [ (iu? t )(iu p rt ,)D a P(u rt -u rt>: t-t') (70) 



z Jrtt 1 

Here the symbol (F[u\) means (F[u + Su])s u and we 
have used the symmetry A"'' = A@ a . We denote 
A al3 '"f{u) = 9 7 A Q/3 (u) etc.. This can be rewritten as: with: 



E^(urt - u rt >,t - = Rll r A(A a ^\u r t - u rt ')) - ff«» ^ R% rt „{A a ^ s (u rt - u rt »))) (71) 
= -^A^e^ (Urt - u -'*' ) (^.i?o,t-t'-^e-' Ft "- (c ' ' - c 'o, t - t ')^_ 5 (t_ i ') /* K.R , T .Ke- K - {Co '°- Co ^- K ) (72) 

D a0 {urt-u r t>,t~t') = (A Q/3 (u rt -u rf )) = ^A^ 3 e- if - (Co '°- c ^*-*' ) - K e lK(urt - u '"' ) (73) 



One can check that the exact FDT equalities are sat- 
isfied (for t > t')\ 

d v D af \u rt - 1^, t - = -T£ Q/3 (u rt - t - (74) 

where the time derivative acts only on the explicit time 
dependence (i.e second argument) of the functior£3 

In the above effective action to first order we will keep 
only the relevant terms, namely: 

(i) the disorder, by definition: 

A°£ = D af3 (u rt - u rt ', oo) (75) 

(ii) the thermal noise, by definition: 

S(VTU = \ JiD^iO, r) - D»P(p, oo)) (76) 

which will be in practice a diagonal tensor 8{r]T) a p = 
5(r]T)5 a p. 

(iii) the friction coefficient, by definition: 

S Vaf3 = lim zc^E^O, w) = - / r£ Q/3 (0, t) (77) 

Jr>0 

Let us note that in the same way that the disorder 
term D generated in the effective action can be split into 
a time persistent (i.e non local in time) operator (dis- 
order (i)) which become relevant at T g and an operator 
local in time (thermal noise (ii)), the term E gives a local 
operator (friction (iii)) and an additional time persistent 
kinetic term (iv) not written above. In the present ap- 
proach these terms are directly related via the FDT. and 
we do not need to consider (iv) separately, though it may 
be important in other cases. 

Using the above FDT relation, integrating by parts 
and using the symmetry C r ,-t = C r ,t one immediately 
finds that TSrj a /3 = 5(r)T) a p, i.e that temperature is not 
renormalized: 



ST = (78) 

The non-renormalization of temperature holds to all 
orders and is guaranteed by the FDT relations. 

Thus, to first order one has simply to consider the dis- 
order: 

A« = A K e- K -( Co >°- Co <*=^ K (79) 
and the correction to the friction: 

1 f+OO 

<W = iVAf / dt( e -*-«*>.o-<*, t ).* 

K Jo 

_ e -K.{C 0fi -C 0%t=oo ).K\ 

Using that A^f = A K K a Kp, A K = gT 2 for the model 
of interest and isotropy one gets that Sr\ a p = Sr]S a with: 

Sr, = 9 — V K 2 [ +CO dt{e- K -^- c ^- K (80) 
2 K Jo 

_ e -K.(C 0l0 -C , t = ao ).K\ ( 81 ) 

which becomes infrared divergent below T g . We now 
turn to the full result up to second order and the renor- 
malization. 



B. RG equations to second order and calculation of 
the dynamical exponent 

The calculation of the effective action T to second order 
in Si n t is performed in the Appendix. We now summa- 
rize all the results up to second order, specializing to the 
case of interest here, i.e model (||) with Ak = gT 2 . We 
have also introduced the dimensionless disorder strength 
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g = go 1 . We have specified a short-scale regularization, with 0(0) = 1. 
which we keep as general as possible (see Appendix for One has: 

all details): 

fly(r,i,o) = 2(Cy(0 > 0,a)-C ii (r,*,o)) (83) 
Q 3 (?,t,a)= T J ( yyV ?f (82) =B iJ (r/a,t/a*,l) (84) 

x \ c i\ Pij e ° liq M * + c 6e^ e Cm<1 M ') With these definitions one finds: 



9R = e 



-i*-B(0,t=oo,«).JC a -2 ( - + ^ fl -2 / ( e Jf.fl(r,0,a).Jf _ 2e -^.B(r,0,a).if ^ 



2\ -if.B(r,0,a).K 



(85b) 



A 66 ,h = A 66 + \fa-*T-£jj^K^ - \{K.r) 

An.fi = A n + -fi^T^ f {\k 2 t 2 + \{K.r) 2 )e- K - B ^- K (85c) 



K 

i r+oo 

r) R = r) + -Tga- 2 V K 2 / dte -^(°>M).K (85d) 



In the first line K' ^ K. The last line is our previ- 
ous first order result (pp|). Studying the variation with 
respect to an infinitesimal change of cutoff a — > a' = e a 
leads to the following RG equations, derived in detail in 
the Appendix: 



dA 66 (l) 
dl 

dl 

dl 

Mi) 
di 



A m {4>)T~g{l) 2 (86a) 

Axri^mi? (86b) 

(2- Kl K 2 T)g(l)+g(l) 2 A g (0) (86c) 

A v (4>)Tg(l) (86d) 



with T c — 2/kiK 2 . The amplitudes Ai{<f>) in general 
depend on the details of the regularisation and are com- 
puted in the Appendix. When evaluated at T = T g how- 
ever, they have a simple form: 



A 66 = ^K 2 e s ^(2I a (a)~h(a)) 

A 11 = ^K 2 e s « c M(2I Q (a)+h(a)) 
A g = 2^ c ^{h{a)~2h{a/2)) 



(87a) 

(87b) 
(87c) 



A„ = ! e T+«(«if2i(* )ii ffc )? = Bv (87d) 
2 c 66 c\i 

One sees that all cutoff dependence drops out of the uni- 
versal ratios determining the critical exponents. 

The dynamical exponent z below T g is determined 
from: 



din 7y(0 
dl 



At the fixed point g(l) — g* = —2r/A g with r = 
(T — T c )/T c one has rj(L) ~ L^~ 2 where the dynamical 
exponent is: 



z-2 = -2r T ^ 



12e 7 r c 66 gas en 

A g 2I (a/2) - I (a) c u c u + c 66 



using K 2 T g = Stt/(c u + Cgg 1 ). This yields our result 
for the dynamical exponent to lowest order in r: 



2-2 



3e 7 



(2 + ")(jrl)^ 
2/ (a/2) - I (a) 



(89) 



with a(T 5 ) = a g = 2(c n - c 66 )/(c n + c 66 ) = 2(1 + 
ex)/ (3 — cr) where a is the Poisson ratio defined as usual 
as a = 1 — 2c 66 /cn = A/(A + 2/i). See fig. 3 for a plot of 
this exponent 9. 

One checks that in the case of the N = 1 Cardy 
Ostlund model one recovers previous result. Performing 
the substitutions explained before equation (^), with 
g* = 7-/(271"), B co = e 7 /c, r g co = 4ttc, one finds: 



z co - 2 = 3% co B co = 2eV 



(90) 



which is the result of Goldschmidt and SchaubBEl] 
(they use yfc = e 7 /2). 

One can obtain further dynamical quantities near the 
transition. Explicitly integrating the RG equation ( |3~l| ) 
one has: 



I' dlg(l) = ±-Hl + X {e 2Tl -l)) 
Jo B g 



(91) 
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This yields the diffusion coefficient D = c/rj at scale 
L = alnl. We denote Dq = D(a) the bare diffusion 
coefficient Dq = D(a). Exactly at T g one finds: 



D{L) 
D 



exp(-BT f dlg{l)) = T. ~ \n(L/a) r e (M) 

J0 ^ + ln(Zo/a)) 



with £o = aex p(l/(.9o^g)) is the characteristic length 
at T g and the exponent = lim r ^ + { z — 2)/t is given by 
the formula (|39|) above. 



Above and near T g one has: 



D g B g 



D eXp( -^ ln(1 ' 2\r 



Finally below T g one has: 



- = d 
D { 



9oB g , 2r i 



2t 



I)) 6 



)) 



,(1-z)l 



(93) 



(94) 



One can also study the velocity v of the system in rp^. 
sponse to a small force /. We folloW|-t|he analysis ofEil 
based on stopped RG arguments as ir£3. Above T g one 
has v ~ fif with fi = D (Einstein relation), and D was 
obtained in (p3|). In the glass phase T <T g one iterates 
the RG until a scale L ~ \fc[f and then match the RG 
to an asymptotic ohmic regime. One thus obtains the 
V — f characteristics as v = fi(L = yjc/f)f. For r > 
one has: 



Mo/ 



(95) 



with /o = c/o? and / c is the effective critical force. 
Since we work at T > there cannot be a true criti- 
cal force, but one can always define an effective critical 
forceEI. In the case of vortices in superconductors it is 
gives the critical current j c . It is also possible to relate 
f c to the Larkin length R c — a(2r / 'goBg) 1 ^ 2 ^ which is 
similar to R a in this single cosine model. The relation is 
fc/fo = (a/Rc) 2 - Thus for weak disorder one has: 



Mo/ 



Exactly at the transition r = one finds: 

/ 



(1 



)) 6 



(96) 



(97) 



with f c = c£ 2 . 

Note that we have used near-equilibrium arguments 
which assume the irrelevance of KPZ type terms near 
the transition and ignores possible violations of Einstein 
relation. A rp«re detailed analysis goes beyond the scope 
of this papeiH. 



V. CONCLUSION 

In this paper we have studied the problem of a triangu- 
lar elastic lattice on a disordered substrate, excluding dis- 
locations. This problem is interesting in relation to vor- 
tex lattices in superconductors, friction of surfaces, mag- 
netic bubbles. We have constructed the N = 2 compo- 
nents model necessary to describe correctly the triangular 
lattice and we have studied both the statics and the dy- 
namics using several methods of renormalization. These 
methods have yielded consistent results. We have stud- 
ied several regularizations and showed explicitely that the 
amplitude ratios which determine the exponents become 
independent of the rcgularization procedure at T g . 

We have obtained that there is a glass phase for T <T g 
where disorder is perturbatively relevant and yields to 
behaviour qualitatively similar, but quantitatively more 
complex than the N = 1 component version of this model 
studied so far. We found that the glass phase is described 
by a a plane of fixed points, parametrized by temperature 
and the Poisson ratio. We obtained that a u 2 ~ A\ In 2 r 
growth of the static displacements and computed the am- 
plitude A to lowest order in T — T g . We found that it 
is a universal function of the Poisson ratio. We showed 
that the asymptotic behaviour of the correlation function 
in the glass phase is isotropic ut ~ Ul, but also found 
a universal subdominant anisotropy — u 2 L ~ A2\nr. 
These behaviours are reminiscent of the behaviour for the 
Bragg glass in d — 4 — e or in d = 3 obtained from the 
GVM, with the difference that u 2 ~ A^lnr in that case 
and !iy — u L ~ est. 

We have also studied the equilibrium dynamics of the 
model. We have obtained the dynamical exponent z, 
which is also a function of temperature and the Poisson 
ratio. This value is larger than the one for the N = 1 
component model indicating that the system is more 
glassy. We also computed the effective critical current 
in the glass phase and related it to the Larkin length. 

Finally, we also reexamined the statics of the N = 1 
Cardy-Ostlund model (planar random field XY model). 
We obtained an expression for the amplitude of the 
u 2 ^ A In 2 r growth of the displacements which seems 
compatible in order of magnitude with the numerical 
simulations, though more extensive simulations would 
be useful. Another prediction is that the distribution of 
rescaled displacements u/ lnr is a gaussian at large scale. 
We propose that this could provide a new and non trivial 
numerical check of the validity of the replica symmetric 
Cardy-Ostlund RG. 

We thank T. Giamarchi, L. Cugliandolo, J. Kurchan 
for useful discussions. 

Note added. 

While this manuscript was under completion we re- 



ceived a prepr int by C. Carraro and D.R. Nelson, cond- 
mat/9607184| who study the same model (g). Their re- 
sults arc different from ours. Their static RG calculation 
does not take into account angular integrals and thus 
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shows no dependence on the Poisson ratio. The two cal- 
culations may coincide only for the case a = 0, but a 
direct comparison was difficult since they did not com- 
pute the In 2 r amplitude. Our result for z, obtained by 
carefully using the same regularization in the statics and 
the dynamics, also disagrees with their result. 



APPENDIX A: REGULARISATION AND 
UNIVERSALITY 

In this appendix we detail the method we used to 
renormalise eqs. (|8q), with a general regularisation func- 
tion </>. This is interesting for several reasons: first as 
we have seen it is crucial to use the same cut-off pro- 
cedure in the dynamics and static RG approach. Sec- 
ond, using a general cut-off gives a well-defined method 
to determine the universal quantities, and moreover it 
might be crucial to avoid the usual logarithmic approxi- 
mation when comparing results with experiments (or nu- 
merical calculations) : for the case of He film, see the 
discussion of Nozieres and GalletH. The procedure we 
develop originates from the momentum renoxmalisation 
of the XY model by Knops and Den OudenO. We use 
here a method that allows a real-space renormalisation : 
thus direct comparison with the hard core procedure is 
possible. 

In order to study a change of the cutoff in eqs. J85| ) 
we introduce: 



hij(r,t, a) 



aT 



(cn 1 



d 2 q <j)'(aq) 
(2tt) 2 q L 

-1 pL -cnq 2 fj,t 



,iq.r 



(Al) 



'66 ± ij 



Note that the above integral for hij(f,t,a) has no in- 
frared divergence in q. First we notice that (with k\ 
defined in (|l4|)) : 



KihijiO, 0,0)^ = -kxTK.K' 



and also: 



hij (r, t, a) = hij (0, 0, a) — -ad/daBij (r, t, a) (A2) 
= hij (0, 0, a) + ]f.V P B i:j (r, 0, a) (A3) 



= h^ (0, 0, a) + td t Bij (0, t, a) 



(A4) 



Then upon a change of the cut-off parameter a — > 5 = 
ae l , we have the following change of cut-off function to 
first order in /: 



4>{aq) = <fi(aq) — qla<fi'(aq) 



(A5) 



This gives: 



Let us focus for ex ampl e on the d efini tion of the renor- 
malised Ag6,_R (eq. 85b ). Using (A6) in it we obtain 



A A i 1 9 2 r r 2l(2~TK 1 K 2 ) f 

A 6 6,ii = A 66 + 2 -I-* e / j \ ' 

K r 



-K.B(r,Q,a).K 



(A7) 



-K 2 r 2 



(K.rf (1-2/ K.h(r,0,a).K) 



which can be put in the same form as the original equa- 
tion (85b) with the following definition of the running 
coupling constants: 



5(0=5e iP-T^) 

dA 66 (l) 



dl 



A 66 (ct>)Tg(iy 



(A8) 
(A9) 



The scaling amplitude (<?!>) depends on the cutoff 
function but not on a which can be eliminated by rescal- 
ing r/d — > r : 



Mf>{4>) 



A" 



(A10) 



x K.h(r,0,l)-Ke- KB ^ K 



Note that this integral over r is by construction finite 
and convergent in the infrared and ultraviolet. 

Similar ly from (85c, B5d) one obtains the equations 
(^6t|j86dD with the RG amplitudes 



A 11 {4>) = -Y,j^-K 2 r 2 + \{K.rf 
x K.h(r,0A)-Ke- KB ^ K 

r + oo 

a^^-Tk 2 dt 

il J ° 



(All) 



(A12) 



Bij(r,t,a) = Bij(r,t,a) + 2l(hij(r,t,a) - hij (0,0, a)) (A6) 



The same method can be used for the RG equation of 
g in the case t = oo: 

Bij (0, cxd, a) — Bij(0, oo, a) + 2l(hij (r, oo, a) — hij(0, 0, a)) (A13) 
= Bij(0, oo, a) + 2lTnx5ij (A14) 

W e us ed hij(r, t = oo, a) = 0. Thus we obtain the RG 
eq. ( |S6c| ) with 

A g (<t>) = 2 f (K.h(r, 0, i), K 'e K - B{r W- K ' (A15) 

+K.h(r,0,l).Ke^ K - B ^- K ) 

As we have seen these RG amplitudes depends in a non 
trivial way on the cut-off procedure (here the fonction <ft) . 
However they enter (for example) the expression of the 2- 
point correlation functions, whose asymptotic behaviour 
should be universal at T c . Thus we have to compute 
the values of the amplitudes Ai(<f>) at T — T c . As we 
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will see they involves the asymptotics of the propagators o. /-n ^ _\ ^ jit 
B«(f,0,a) (|f| »a) and%(0,f,a) (i>a 2 ). 

First we focus on the r-propagator, which gives also the 
asymptotic interaction of the Coulomg-Gas (see (|l4|)): 

V(f, a) = Bij (r, 0, a) = E{r)5 l3 + F{r)hfj 



Ki I In 



4t/x 



- 7 - 4ttC*(0) (A25) 



-4- ( c ii lo g c ii + cee 1 iogcee) 



where -F(r) is given by 

F(r) = 2T±-(c^ - c^) J dq^-Uqr) (A16) 

The other coefficient is given by 2E(r) = Tr.B(r) — 
F(r) where 



Tr.B{r) =^B ii (f,0,( 



(A17) 



?(c n 1 + c 66 1 ) / ^(rVa)ln 



The large r behaviour of the integral follows from the 
formula (|r| ^> a): 



1 



9 z 



(A19) 



where C{4>) = <f~ J (Pf^r') ln(r'). The asymptotics of 
the interaction is then : 

E(r) w 2Tki [ln(r/a) - 2ttC*(0)] + Tk 2 (A20) 
D(r) w -2Tk 2 (A21) 

Now we turn back to the time propagator, which is 
given by 

Bij(0,t,a) =TSij [c^F(cu,t,a) +CQQF(c 66 ,t,a)] 

(A22) 

where we define the function 
d 2 q 4>(aq) 



?>■'■">= / (^yM l-^ (A23) 



47T 



d 2 r(j>{r)Ei 



Act 



(A24) 



where -Ei means the Exponential-Integral function which 
verifies (7 is the Euler constant) : 



Ei(-x)= 7 + Hx) + ^2^ 



n=l 



Taking the limit of large time in ( A22] ) using the above 
equation, we find 



We are now able to find the expression of the RG am- 
plitudes Ai(4>) at T c . We will give a detailed calculation 
for A „((j>), the others being similar. Using the property 
(|A|) : K.h(Q,t,l).K = -2T/T C + td t (K.B(0,t,l).K), 
we can write A r) (4>) (eq. (A12)) (after an integration by 
part) as 

M4>) = -\ E (2(1 - 1) [ dte-iWA" (A26) 

-2 U e -iK.B(0,t,l).K 

J 

This expression is valid for T > T c ,&s well as for T = T c 
provided the first term is se t to . With the asymptotic 
form of the time correlator ( jA25| ) : 



e 2 



|K.B(0,i,l).K 



ex P [--TK 2 (- Kl (7 + 47rC(0)) 



with 
T K 

+ ^~( c ii logcn + logc 66 ) 
we obtain the expression at T c : 

A n = \^2K 2 fi- 1 T K (T = T C ) 

K 



(A27) 



(A28) 



We now t urn back to the other amplitudes 
( All , A10 A15 ). By the same method, and using the 
asymptotics (A21) : 



e -K.B.K _ r -4£ e 2a(K.r) 2 e -a+$TrC(<j>) 



(A29) 



where a has been defined in section II : a = k 2 K 2 T 
we obtain at T = T c : 



-466(0) 



1 



i+8ttC(0) 



E 



2 7; 



K 2 d0 



(A30) 

_I cos 2(^ e 2«cos 2 W 



An{<f>) = = ; e- a + a "CM 



E 



K 2 d9 



1 + l cos 2(^ e 2acos 2 W 



For A g the integral split into two terms and we need 
another asy mpto tic expression (the first one follows di- 
rectly from QA29D): 
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,\K.B.K' _ r -2^ e -f +4*C{<l>) e -2 a {k.r){k'.?) 



Thus 

A„ = 2e ~f + 4 - c W 



(A31) 



271 



^0 I j_ e -2acos(e)cos(e+27r/3) _ g-acos^) 



Performing the angular integration one finds the for- 
mulae (p7j) in the text. 



APPENDIX B: STATIC RG CALCULATION 
FROM THE EFFECTIVE ACTION 

In this Appendix we obtain the perturbation theory 
result (^5|) from the statics replica calculation, using the 
effective action method. 

The replicated partition function is Z n = Tr n e~ s with 
S[u] = So[u] + Si n t[u] and: 



ab 



ab 

We will study the general model: 
R(u) = J29Ke iK - u 

K 



(Bl) 



(B2) 



In the case of the triangular lattice, the sum is over 
the six reciprocal lattice vectors K = 1, 2, 3, —1, —2, —3. 
With no loss of generality R{u) is an even function and 
thus qk — 9—K and gx is real. Note here the re- 
dundancy in the sum under the symmetry (K «-> —K, 
(a,b) <- (M). 

The original model (H) we want to study is: 



R(u) = g cos(if.i 



(B3) 



if=l,2,3 



where we denote by 1,2,3 three reciprocal lattice vec- 
tors. This model is obtained as a particular case of (B2) 
for g K = g/2. 

To obtain the effective action T[u] to second order in 
disorder, one has to study all possible Wick contractions 
of the two vertex operators, with all one particle reducible 
diagrams being absent (see below): 



~ f R « ~ u r + S < ~ Su b r )R(u c r , - uf, + Su c r , - 

Z J rr f 

abed 

(B4) 



2r< 2 > - (R( Zl )R(z 2 )) - (R( Zl )) (R(z 2 )) (B5) 
-(R'(z 1 ))(z 1 z 2 )(R'(z 2 )) 



In Fourier it reads: 



R{q 1 )R{q 2 )e-^^ ) e^^^ ) 



(B6) 



The result is: 



Jrr ' abcdKK' 



) e iK'.{u" T ,-u d r ,) 



(B7) 



with 

-wrab.cd/ \ 

V KK> ( r ) = 9K9K 



x ^ e -K.(C c r a + Ct b -C c r b -C$ a ).K' _ ^ 

+ K. {C c r a + Cf - Cf - C d r a ).K') (B9) 

The points r and r' are constrained to be within a small 
distance. The contraction of the points r and r' result 
in new local operators which will correct the effective ac- 
tion. The coefficients can be computed using a standard 
short distance expansion. The relevant new terms are of 
two types: 

(i) the fusion terms: 

they are of the form: 



r? 



(BIO) 



ib K 



These fusion terms can be split in two types, 5gx = 
&i9k + S 2 gx, the same K vector Sigx and different K 
vector 5 2 gx- 

(il) same K fusion terms: 

The terms ( K' = K, c = b or d = a) and (K' = -K 
and c = a or d — b) lead to: 

hS* = 11^ ( V K C K b (r) + V^(r) (Bll) 

(there are of course four other terms obtained by the 
symmetry (K <-> —K, (a,b) «-> (b,a)) which correct the 
corresponding symmetric term). 

This correction can then be computed explicitly. One 
finds that for a ^ b performing a sum over c ^ a,b: 

S l9K = 2{n-i)e-i KB "^ K {g K f ! c -i KBI 'l K Air) 
**> ' (B12) 

with 



The contractions are of the form: 



A(r) = 1 - e -K-C<r)-K(i + K.C{r).K) (B13) 
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Th is y ields, upon setting gx = <?/2, part of the equa- 
tion df). 

(i2) different K fusion terms: 

These terms are a peculiarity of the triangular lattice. 
This other correction comes from a = b and c = d, or 
a = d and c = b. 



S29K+K' = ^ / (V^(r) + V«%,{r)) (B14) 



For a given K3, for each (K\, K2) such that K\ +K2 = 
K 3 one finds: 

S 29 k 3 = 9K 1 9K,e-*"*- B ' h W / ^-"W* (B15) 



This yields, upon setting = $/2, the remaining part 
of the equation (85), taking into account that two couples 
{K\,K?t) and (K 2 ,Ki) correct gK 3 - 

(ii) annihilation terms: 

Finally the subcase K' = —K gives a correction to the 
quadratic hamiltonian. It gives the local operator: 



r > = (i)/EBi 

1 Jr ab K Z 



(B16) 



x [(r.V)K " u b ).Kf{V«%{r) + V££>)) 



ib.ba 



which gives: 

T 2 = -d a u a 1 dpu b s ^{9KfK 1 K 6 [ r a 



773 e 



-K.B(r).K 



(B17) 



This should correct 
1 - 



— 9 a «"9|S«j((Au - ^m)SarfSl3S + A 66 ((5 ct(3 5 7( 5)) 



2T 



Let us define a tensor: 



(B18) 



Ta-y,l3S — ASapSyS + — (SayS/lS + S a gS/3^) (B19) 

Then one finds: 

4 _ 3 1 1 1 

o 4 4 Z 

(B20) 

Using this algebra it yields, upon setting gx = 5/2, 
the corresponding part of the equation (Bq). 



APPENDIX III: DYNAMICAL EFFECTIVE ACTION TO SECOND ORDER 

In this appendix we derive the perturbative expression of the effective dynamical action to second order in disorder. 
We then identify the relevant terms which correct the bare disorder and lead to divergences below T g . This procedure 
amounts to a short distance, short time operator product expansion. Note that the operators are local in r but non 
local in time, which makes the expansion more involved. .— . 

The effective action to second order in the interaction term isE2l: 

-2T^[U] = (S int [U + 6Uf) su - (S int [U + 5U]f su - ( ^[^+ 6U] ^^ 68^^ + SU] ^ (1) 

with U — (u, it) and SU = (<5it, 5u) and a gaussian average over SU is performed using the bare quadratic action 
So + S2 in (p8|). The last term merely ensures that all one particle reducible diagrams be absent. 
One thus has to study all possible Wick contractions of the two vertex operators: 

iU^iKl^Aaifoiu^Jx ~ U rilti )A Q2/ 3 2 (u r2t2 - U r ^t' 2 )iUrX iil rtt' 2 ( 2 ) 

imposing that at least two contractions joining the two vertices 1 and 2. 
A tedious calculation gives, for the iuiu term (in short notations): 

-2r = 2( i < 1 )^ 2 2 )(A ai/3l K 1 - u t[ )A a2 p 2 . nS (u t2 - ^ 2 )) c <4(<; - RlP t \) 

+ i(z^)(z^)(A Ql/3l ( Mtl - u t ,)A a2f32 .„ 5 (u t2 - u t . 2 )){RZ\ ~ Rl£)(RZ\ ~ R'ilO + 

( lU ^)(zu%)(A aiPl .Au tl - u t[ )A a202 . n (u t2 - u t ,))(Rl% - R^(R S t X ~ 4%) 

+(A ai ^ 5 (u tl - u t ,)A a2 p 2 . n {u t2 - Ut -)) c ( i ^)K 1 )< t2 (<? 2 - R^ t 2 2 ) (3) 

Using the assumption of time and space translational invariance it can be put under the form T = 
-±(iu? ti )(iu p r+r , t2 )5Af. One finds four terms: SA = £4=1,4 8A® f : 
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SA { r ] ] = 2R sx {T 2 ,Q)W p {T ll r'){A px , lS {u r+r , M - u r+r , M - T2 ) 

\^Oip{^rt\ ^r,t2— Ti) ^api^rti ^r.t2 — T\ — T2 )])c 

SA^} = ^5(r")R-'P(T,~r r )R 5x (T / ,^r')(A a0 , 1 s(u r t 1 - U r t 2 ) 

[ApX^r+r'jx-r ~ Ur+r' ,ti—T') + ^p\{ u r+r' ,t 2 -T ~ Ur+r\t 2 -r') ~ 
A p A(Wr+r' : ti-T — Ur+r' ,t 2 —r') ~ ^px{u r +r' ,t 2 —r ~ u r+r',ti-T')J) 

5A^ = R lp {T2,r')R sx {Ti, -r')[(A a p-s(u rttl - u r A 2 -r 2 )Af}x n (u r+r ^ t2 - u r + r ' m-ti)) ~ 

( A ap; S {U rt t t — U r J 1 -Ti-T 2 )Af}\^ 1 {u r + r iJ 2 - l* r +T J ,tl-Tl)) — 

A a p-s {Urfa — u r j 2 - T2 )Ap\;y(u r+r ij 2 — U r +r' ,t 2 -ri-T 2 ))\ 
5£$ = S(r")R^(T 2 ,O)R sx (T 1 ,-r')(A a0lS (u rM - u rM ) 

\A\p-~j (ii r _)_ r '^ 1 _ Tl u r j rr ' f l — Tl — T2 ) A\p-j (ii r _)_ r '^ 2 _ T1 u r ^- r ' i 2 _ Tl _ T2 )] ) c (4) 

We have written these terms in that form for simplicity, but one must keep in mind that in addition they must be 
symmetrized under a — > [3 and r — > —r when necessary. Note that, to avoid making the notations even more heavy, 
we have deliberately chosen in this Appendix not to write systematically all space and time integrals and we assume 
that the reader can restitute the necessary integrations. 

We now identify the new operators generated in a short distance expansion. As discussed in the text we will not 
pay special attention here to the kinetic term (i.e the friction to second order as well as the contributions to the time 
persistent kinetic term) but rather focus on the contributions to the disorder. Note also that it is also possible to 
use a supersymmetric perturbation theory formalism to study the. dynamics of this model. However since we may 
also be interested in the future to cases which are not potentialEj and to which supersymmetry does not apply as 
conveniently, we will use a more pedestrian (and painful) approach. 

The above terms (Q) give either a contribution of fusion that we write under the form: 

SA a ^e iK ^- Urtl ~ Urt2 ' e~^ K - B °< tl - t2 ' K (5) 

r,t±,t 2 

( We have defined = 2(C^ - C«£), identical to definition (||) and we use that ti-t 3 +n = B o ti-t 2 + 

Ofa/fa-ta))), 

or, they give a contribution annihilation of the form: — ^(iu" tl )\7 2 (iii^ +r , t2 ). This analysis parallels the one of the 
preceding section from a static replica calculation. 

(i) different K fusion terms 

The terms which give (from the Coulomb gas analogy) the fusion of the type K\ + K 2 — ► if 3 are the last two terms 
of <5A( 2 ) and the first term of SA^ 3 \ The last two terms of SA^ give: 

8*k+k> = ^ 7 ^A^A^,ii^ >T ii^ iT ,e x ( 2Co -°- c --- c --') K ' (6) 

where a factor KBQ, tl -t 2 K' has been absorbed in front to make exp(— 1/2 (if + K 1 ) 2 £?o,ii-t 2 ) appear in front, as 
required by the definition (Q) of A^f +K , , and 



SAf_ K , = iif 7 if 6 A^A^ A ,^ T iii A nT ,e- K -( 2Co . - c ^--- c -.--')- K ' (7) 
Adding to these terms the first term of <5A' 3 ^ gives a total contribution: 



SAf +K , = ^K 5 A a K (A p K x , + A^,)i?l P r , r i?i A r , r 'e K -( 2Co ' ~ c — + 

^ifaA^A^^WH^e^^ 00 - -^--^.--')^' (8) 

Using now that A^f = Ai{K a Kp, the symmetries A°^f = A^f , the r — > — r symmetry of R and C, and the r — > — r 
symmetry of C, one finds that all three terms combine into: 

SAf +K , = A K A K ,(K a K + K a K' )(K.R r>T .K')(K.R rtT ,.K')ei K -^+ B ^') K ' (9) 

Putting together the contribution of (K,K') and (if', if) one has finally: 
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SAf +K , = A K A K ,(K + K') a (K + K') [ {K.R r , r .K'){K.R rtT ,.K')e^ B ^ +B ^ K ' 

J r,r,r' 



(10) 



This expression can be simplified using the FDT theorem (^4|) 9tC™f = Ti?"f . Integrating by part and noting that 
B r ,r=o = B stat {r) one finds: 

6Af +K , = A K A Kr ^(K + K%(K + K%Je K - B ^ K ' (11) 

we have also assumed K.B r T=00 .K' — — oo. This is yields the result (85a) if one uses Ak = gT 2 . It coincides also 
with the result of the statics (see previous section). 
(ii) annihilation terms 

The annihilation term (i.e the first term of SA^ for K' = K and the last two terms of 5A^ 2 ^ for K' = —K ) gives 
(assuming the K — » —K symmetry): 



e -K.{2C 0<0 -C r , _r-C T , _ T ,).K 



6Af(r) = W K 7 K s (6(r) AfA" K x f Rl p rl>T R%, r , 

A% p A K x R^ T Ri\ T ,e- K ^ 2Co ^- c ^~ c ^ '1- K (12) 
Thus assuming all symmetries it yields: 

(13) 



is Jt.t' 



-K.(2C Q , Q -C r , T -C T _ T ,).K 

K 

The time integrals can again be performed using FDT. One finds: 



/ 7 V) = ^(AK) 2 e- K - B -° * (W) 



j, 2 

Since, by definition 5A2°(q) ~ — T[(An — A 6 6)g 7 9<5 + A 6 6(? 2 5 7 5], one gets: 



if 



{SA n - SA 66 )^(6 ia 5s f 3 + 6 lfJ 5 Sa ) + 5A 66 5 af3 § lS = \ I r a r» f\r) (15) 



which yields equations (85b 85c) when using Ax = gT 2 . 
(Hi) same K fusion terms 

Finally the same K fusion terms are the most delicate. There are a priori jour distinct contributions. Let us start 
with the first one, which as we show vanishes identically: the last two terms of 8A^ give, respectively: 

8Af = K s A^J2K A K'^^Ari e ~^' Bo,T1+T '' K ' e ~ K ' lCr ' -^- Cr ^- K ' (16) 

K' 

and 

8A°f K = K^A^Y^K'sA^R^R^e-^'- 3 *^^ (17) 

K' 

However these last two terms give zero identically. Indeed using all symmetries and definitions and remembering 
the necessary symmetrization one finds first that the last two terms of SA^ give: 

SAf = A K A KI (K , .R r , n .K){K'.R r>n .K)e-^ K '-^+^ K ' {K a K' & - K K' a ) e - K -^- c ^- K ' (18) 

K' 

Each term vanish because of the symmetry t\ — ► r 2 which makes the summmand over K' odd under K' — > —K'. 
The three remaining contributions are: 
(i) The 6A^ term which gives: 
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6Af = -2A^^^;^A^i?^ r2 ^ Ti e-^'- s o-2^'( e ^(^, T1 -c r , T1 _ T2 ).K' _ ^.(c^^-c^.k^ 



K' 



(ii) the first two terms of SA^ which give a contribution: 

5Af = -K^KsAf ^ A p K x ,Rl p rT R s _\. T ,e-i K '- B ^'-r- K ' coshK.(C r ,. T - C T ,- T ,).K' (19) 

K' 

(iii) the term SA^ gives: 

SAf = KsA^YsK^k'^R 5 ^- 1 **'-* ^ - e^cv.-i-Cr,-^).*') 

if 

Using FDT we will show that the last two contributions will almost cancel. First we need to rearrange the first 
one. The SA^ term gives: 

6A$ = -J2(K a K^ + K' a K p )A K A K ,(K\R a , T2 .K')(K\R r ^.K)e-^ K ' B ^ K ' 

K' 

( e K.(C r , T1 -C r , T1 -T 2 ).K' _ e K.(C r , T1 + T2 -C r , T1 ).K'j 

We now split this term in two parts: 

(JA^ 1 ) = -YS.K a K' p + K' a K ) A K A K , [ (K'.R , n .K')(K'.R rtTl .K) 

Xi Jt 2 >t 1 
e -\K'.B 0tT2 .K'+K.(C r , T1 -C r , T2 - T1 ).K' 

which can be changed into, using FDT: 

5A^ 1] =-J2(K a K' + K' a Kp)A K A K , f (K> ' .R r , T2 _ Tl .K)(K> \R r , Tl .K) 

W Jt 2 >Ti 



j^i Jt 2 >ti 
±K'.B , T2 .K'+K.(C r , T1 -C r , T2 - T1 ).K' 



e 2 

+ (h^ K » K P + K '» K ^^A K , f (A". J R r>T1 .A')[ e -^'- B o.r a -«"+^(Cr, T1 -C P . T2 _ T1 ).Jf' ] » = 



(20) 



The first integral can be rewritten using r = n and r' = r 2 — n and becomes symmetric in r — > r'. It is still odd 
in if' however and thus is identically zero. Thus: 

SA^'V = (hj2(K a K'p + K' a K )A K A K , f (^.^.^[e-^'-Bo^^'+^-^^-c,,,,^).*:']^ (21) 

The second piece is: 

5A%' 2) =-Y J {K a K' + K' a K )A K A KI {K'.R o ^ 2 .K'){K'.R rtTl .K)e-^ 

K' 

( e K.(C r , T1 -C r , T1 - T2 ).K' e ( Ti _ ^ _ e K.(C r , T1 + T2 -C r , T1 ).K'j 

This last term can be transformed, changing t\ — > Ti + t 2 in the first integral, into: 

<5A^' 2) = - J^a^is + KK(})A K A K ,(K'.Ro, T2 .K')(K'.R r , Tl+T2 .K - K'.R r , Tl .K) 

K' 

e -\K' .B , T2 .K' e K.(C r , T1+T2 -C r , T1 ).K' ( 22 j 

which can be integrated using FDT into: 



u' Jt 2 



-C r ,T 1 )-K'-iT 1 = oo 
\t 1= 

K' 
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Similarly the first two terms of SA^ give a total contribution: 

6A% } = -^^^A^AKK^-^r.^O^-^r'-^Oe-^'-^-'^^'e^^*--- *--^^' (23) 

K' 

The term SA^ gives: 

5Af =2if Q ^^A K A K^ (^^i? OIT2 .i^0(^•i?r,r 1 •^0e~^'■ Bo ' T2 ■ K 'e" X • (C '''• T ^ C ''• T1+T2) ■ K ' 

K' 

Using FDT this term and integrating by part over ti this term gives: 

8A% ] = -2K a K^ AKA^(if'.i? r!Tl+T2 .X)(/f.i? r!Tl .X')e"^'- Bo -2^' e -^(c., T1 -c., T1+ . 2 ).x' _ 

K' 

'.(C riT1 -C riT1 + T2 ).K'^T 2 = 00 



(-^)2K a K ^A K A K , f (K.R rtn .K')[e-* K '- B °™- K '- K - 1 



\t 2 =0 



K 

and thus: 



SA { x ] = -5A^ + {-)2K a K fj ^A K A K , f (^. firTi . j ft:')[ e -K'-s , T2 .K'-x.(c t ., T1 -c t ., T1+T2 ).if']r2=oo 



(24) 



where r = n + t 2 and r' = n has been used, as well as the symmetry in r — > r' of <5A( 2 ) which allows to write it 
as twice the integral over r > r' . 
Thus, one has: 



6A$ + SAP + 6A$ = (i) + K' a Kp)A K A K , 

K> 

( / (^^i^ r;Tl .^)[ e -|^^Bo,. 2 .^'+^f.(C' t ., T1 -c r , T2 _ T1 ).^f']r2=oo + 

J Ti 

/ (i^^i?o : r 2 .^Oe _ ^'• Bo '^•^'[ e ^•(C., T1+T2 -C r , T J.K']rl=oo) + 

(i)2^ Q ^2 A ^' / (^•i2r,Ti.«")[e-^'- Bo ^- Jf '- Jf 



(C r>T1 -C r , T1+T2 ).-ff'-|T 2 =oc 
Jt 2 =0 



(25) 



K 

This can be rewritten as 



6A$> + 5A%> + 5Af = d)J2(K a K>p + K' a K p )A K A K , 

K> 

J{K\R r ^.K)e- K ' c ^ K ' +KC ^ K ' + ^^ 

- J[K'.R TiT .K + K\Ro, T .K')e- K '-( Co -°- c °^- K ' +K -( c ^- c ^- K ' 

+ d)K a K p £ ^kA k , f (K'.R r , T .K){e- K '- c ^- K '- K '- c ^- K - 1) (26) 

Using the symmetry K' — > —if' one see that the second term vanishes. Using FDT and this symmetry one finds 
that the first and third term cancel. Finally one finds: 

SA^ + 5A$ + 5A$ = -(^)K a K p ]T A K A K , [ e -* W -k> .c r . .K (27) 

The divergent integrals at T g correspond to K' = —K, with the final result: 

5A$ + 5A$ + 5A<$ = -( A )K a K J2(*Kf I e~ K ^- c ^- K (28) 

K> Jr - r 

Note that we have used that C r;T=00 = 0. 
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APPENDIX IV: STATISTICAL SYMMETRY AND 
CORRELATION FUNCTION 

In this Appendix we detail some useful consequences 
of the statistical tilt symmetry. It allows for instance to 
obtain a useful form for the real space Callan-Symanzik 
equations referred to in the text. We consider for sim- 
plicity the N = 1 planar random field XY model defined 
in its replicated version as (52), the generalisation to (|^) 
being straightforward. 

First we derive the Ward Identities of the statistical 
tilt symmetry on the replicated model. Under a local tilt 
u a {x) — > u a (x) + e(x) the correlation generating func- 
tional ZW[J] remains unchanged. This invariance can 
be written 

/ d 2 x < X ) ^ E ^ + E '•(*)) zin) w = 

Introducing the generating functional of the connected, 
correlation W^ n ' [J] and using a Legendre transform^ 
we arrive at the Ward identities for the 1PI replicated 
functional \u] : 



c 
T 



EW( 3 



E 



<jr< n >[u] 

5u a (x) 



(1) 



^From this one deduces useful consequences, e.g that c is 
not renormalized to all order. 

Another useful consequence of this symmetry is to re- 
late the correlation functions of the original model to the 
one of the model with A = 0, as noted ir£3. This is easier 
to demonstrate on the unrcplicatcd model (analogous to 
(|)). Consider: 



G h (x) = ( e *f h (xMx)} 



(2) 



Using the local shift Vu — » Vu — —a where a is the 
local random stress for the N = 1 version of model (|^) 
(one can assume a q to be purely longitudinal), we can 
(up to a redefinition of the random field <fi which does 
not modify the average ) factorize out the A dependancc 
of the correlation function (0) : 



G h (x) 



' 2 la 



rh a h- 



G h (x)\ 



A=0 



(3) 



lnZ(A,T) = 



p 2 TA f d 2 q 1 - cos(g.r) 



(2n? 



Using formula ( A18) we find 



lnZ(A,T) 



p 2 TA 
' 2ttc 2 



In 



(D 



(6) 



(7) 



This can now ben-used to study the real space Callan- 
Symanzik equationED. The two point function behaves 
under a rescaling as: 



G(r, A(0), g(0)) = e~ li *G{re- 1 , A(l),g(l)) 



(8) 



where d* is the anomalous dimension of the operator 
e ipu(x) ( wn j cn was computed e.g in Appendix B o 0). 
The invariance of this correlation function under a change 
of cut-off parameter a — * ae l gives the Callan-Symanzik 
equation: 



d 



dA d 
~dT~dA 



dg d 



d\nr p dl dA dl dg 
Using the above relation into (fl) we obtain 



G(r, a) = (9) 



d 



d In i 



P 2 T /r\ dA(l) dg 3 







Thus we obtain that if the operator e lputyX ^ is multiplica- 
tively renormalisable to all order, the dominant term in 
the large r limit will be to all order an exp(A. In 2 (r/a)) 
term. To first order, using ( |3p| ) we find on the low tem- 
perature fixed line (dg/dl — 0), near T g = Ai:c,g = t/2t: 



In G(r, a) 



- P 2 t 2 In 2 (-J + OQn(r/a)) 



Note that this shows that the dominant term (r 3> a) 
has a gaussian distribution. Thus one obtains permuting 
the limit p — > and the disorder average) : 



B(r) 



-2^r In G(r, a) = 2r 2 ln 2 (r/a) + C(ln(r/a)) (10) 
p 1 



which is the result obtained in Section ( [II B 2 ) by an- 
other method. 



If one consider for instance the regularized two point 
correlation function of the operator e lpu ^: 



G{x,a) = ( e wu(x) e ~*pu(o)} 



(4) 



obtained from the above with for h(y) — pS(y — x) 
pS(y), one has: 



G(r,o) = Z(A,T).G(r,o)| A=0 



(•5) 



where G(r, a) 



A=o 



is the correlation function with A 



set to zero in the model (), and the factor 
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